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Abstract 

In this note it is proposed a class of non-stationary de Sitter, rotating and 
non-rotating, solutions of Einstein's held equations with a cosmological term of 
variable function A*(u). It is found that the space-time of the rotating non- 
stationary de Sitter model is an algebraically special in the Petrov classification 
of gravitational field with a null vector, which is geodesic, shear free, expanding 
as well as non-zero twist. However, that of the non-rotating non-stationary model 
is conformally flat with non-empty space. 

Keywords: Non-stationary de Sitter; rotating and non-rotating cosmological mod- 
els; Kerr-Schild ansatz. 

It is well known that the original de Sitter cosmological model is conformally flat 
Cabcd = space-time with constant curvature R abcd = (A* /3)(g ac g bd - g a d9bc) [1]- It 
also describes the non-rotating and stationary solution. Therefore, the non-rotating 
stationary de Sitter model is a solution of Einstein's field equations for an empty space 
with constant curvature, whereas the rotating stationary de Sitter model proposed in 
Ref. 2 is a solution for non-empty space with non-constant curvature. Because of the 
stationary and non-rotating properties of the original de Sitter space, the non-rotating 
Schwarzschild black hole with constant mass can embed to produce Schwarzschild-de 
Sitter cosmological black hole with two event horizons - one for black hole and other for 
cosmological [3]. Similarly, the rotating stationary de Sitter cosmological universe [2] 
can conveniently embed into the rotating stationary Kerr-Newman solution to produce 
rotating Kerr-Newman-de Sitter cosmological black hole with constant cosmological 
term. This Kerr-Newman-de Sitter black hole metric can be expressed in terms of Kerr- 
Schild ansatz with different backgrounds as g™ dS = g^ b +2Q(r, 0)£ a £ b where Q(r, 9) = 
-(rm - e 2 /2)R~ 2 , and <?*NdS = gf b N + 2H(r,6)£ a £ b with H(r,6) = -(A*r 4 /6) iT 2 . 
Here g^f is the rotating stationary de Sitter metric and the vector £ a is a geodesic, shear 
free, expanding as well as non-zero twist, and one of the repeated principal null vectors 
of g^b , dab an d <7^ NdS > as these space-times are Petrov type D. The expressibility of an 
embedded black hole in different Kerr-Schild ansatze means that, it is always true to 
talk about either Kerr-Newman black hole embedded into the rotating de Sitter space 
as Kerr-Newman-de Sitter or the rotating de Sitter space into Kerr-Newman black 
hole as rotating de Sitter-Kerr-Newman black hole - geometrically both are the same. 
That is, physically one may not be able to predict which space starts first to embed 
into what space. One thing we found from the study of Hawking's radiation of Kerr- 
Newman-de Sitter black hole [4], is that, there is no effect on the cosmological constant 
A* during the evaporation process of electrical radiation. The cosmological constant 
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A* always remains unaffected in Einstein's field equations during Hawking's radiation 
process. That is, unless some external forces apply to remove the cosmological term 
A* from the space-time geometry, it continues to exist along with the electrically 
radiating objects, rotating or non-rotating. This means that it might have started to 
embed from the very early stage of the embedded black hole, and should continue to 
embed forever. It is noted that the Kerr-Newman-de Sitter black hole proposed in 
Ref. 2 is found different from the one obtained by Carter [5] in the terms involving 
cosmological constant. 

The black hole embedded into de Sitter space plays an important role in classical 
general relativity that the cosmological constant is found present in the inflationary 
scenario of the early universe in a stage where the universe is geometrically similar to 
the original de Sitter space [6]. Also embedded black holes can avoid the direct for- 
mation of negative mass naked singularities during Hawking's black hole evaporation 
process [4] . It is also known that the rotating Vaidya-Bonnor black hole with variable 
mass M{u) and charge e(u) is a non-stationary solution. When M(u) and e(u) become 
constants, the rotating Vaidya-Bonnor black hole will reduce to the stationary Kerr- 
Newman black hole. If one wishes to study the physical properties of the gravitational 
field of a complete non-stationary embedded black hole, e.g. rotating non-stationary 
Vaidya-Bonnor-de Sitter (not discussed in this note), one needs to derive a new ro- 
tating non-stationary de Sitter model with a cosmological term of variable function 
A*(u). That is, an observer traveling in a non-stationary space-time must also be 
able to find a non-stationary cosmological de Sitter space to embed, having a similar 
space-time structure with time dependent functions. 

In this view, it is proposed a rotating non-stationary de Sitter solution of Einstein's 
field equations with a cosmological term of variable function A*(u) in this note. Using 
Newman-Penrose formalism [7], a class of rotating metric with a mass function M(u, r) 
has been discussed in Ref. 2, where the mass function is being expressed in terms of 
Wang-Wu function q n (u) [8] as 

+oo 

M(u,r) = Qn(u)r n . 

ra=— oo 

For obtaining a rotating non-stationary de Sitter solution, we choose the Wang-Wu 
function as 

„ (lA _ / A *M/6, when n = 3 

qn[U) ~ \ 0, when n + 3, Uj 

such that 

M(u,r) = Va*(u). (2) 

Then using this mass function in the rotating metric presented in equation (6.4) of 
Ref. 2, we obtain a rotating metric, describing a non-stationary de Sitter model with 
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cosmological term A* (it) in the null coordinates (u,r,6, 4>) as 

ds 2 = |!_!^0j dtt 2 + 2dudr 

+2a ^ ; sin 2 ^ dn # - 2a sin 2 6> dr # 

-R 2 d6 2 - {(r 2 + a 2 ) 2 - A*a 2 sin 2 #} iT 2 sin 2 dcp 2 , (3) 

where R 2 = r 2 + a 2 cos 2 6> and A* = r 2 — r 4 A*(u)/3 + a 2 . Here A* (it) denotes an 
arbitrary non-increasing function of the retarded time coordinate it and a being a 
constant rotational parameter. When one sets the function A*(u) to be a constant, 
the line element (3) will reduce to the rotating stationary de Sitter space-time [2] . The 
complex null vectors for the above metric can be chosen as follows: 

e a = 5l-asm 2 6 5 4 a , 

na = ^S 1 a + S 2 a -^asm 2 eSi, (4) 
m a = -^2R { ~ iasindS l + R 2 ^a+ *(r 2 + a 2 )sin#o 4 }. 

Here £ a , n a are real null vectors and m a is complex with the normalization conditions 
£ a n a = 1 = —m a fh a . By virtue of Einstein's field equations, we calculate the energy- 
momentum tensor describing matter field for the non-stationary space-time as 

T ab = At* 4 4 + 2/?*^( a n b) + 2pm (a fh b) + 2ul {a rh b) + 2u)t {a m b) , (5) 

where 

4 4 

I 1 * = ~ 6KR 2 R 2 { 2rA * {U) + ° Wg A * {U) ' uu } ' P * = KWm k * {u) ' 



P 



■■ 2A -<*V + 2«W*}, (fl-aS)^.),, (6) 



KR 2 R 2 I J ' 6^/2KR 2 R 2 



with the universal constant X = 87rG/c 4 . The quantity p* interprets as a null density 
based on the derivative of A* (it); p* and p are the density and pressure of the non- 
stationary matter field, and u represents the rotational force density determined by 
the rotational parameter a coupling with the derivative of A*(u). That is, for non- 
rotating model a = 0, the rotational density to will vanish, describing a non-rotating 
fluid model for non-stationary de Sitter universe. From (6) we find the equation of 
state (the ratio of the pressure to energy density ) for the non-stationary rotating 
solution as w = p/p = — (r 2 + 2a 2 cos 2 #)r~ 2 for non-zero rotational parameter a. This 
equation of state will take the value w = — 1 at the poles 9 = tt/2 and 6 = 3ir/2, 
showing that the de Sitter solution (3) with non-constant A*(u) describes a rotating 
non-stationary dark energy model possessing negative pressure. 
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The trace of energy momentum tensor T a b (5) is found as 



T = 2(p*-p). (7) 

Here it is observed that p* — p > for rotating non-stationary de Sitter model. The 
energy-momentum tensor (5) satisfies the energy conservation equations T ab . b = 0. 
The verification of these equation may be seen in Appendix A below. The Ricci scalar 
A (= ^g ab R a b), describing matter field by virtue of Einstein's field equations, is found 
as 

A= ^r 2 A*(u)R- 2 . (8) 

Other Ricci scalars are related with p* , p*, p and w (6) as Kp* = 2fa2, Klo = —2 012, 
Kp* = 2 <fin +6 A and Kp = 2 (fin— 6 A. The energy momentum tensor (5) also satisfies 
all the three energy conditions: (i) weak, T ab U a U b > 0, (ii) strong, T ab U a U b > \T and 
(iii) dominant, for a time-like observer with its four-velocity vector U a as shown in 
Ref. 2. It is noted that T ab (5) does not describe a perfect fluid, i.e. for a non-rotating 
perfect fluid, = (p* + p)u a Ub — pg ab , with unit time-like vector u a and trace 

T"( pf ) = p* — 3p, which is different from the one given in (7). 

The existing Weyl scalars, determining gravitational field for the space-time metric 
(3) are obtained as 

r 2 A*(u) 

fa = ^jr^{( r + 2ia cos °) R ~ rR }> 

i a r_ 3 sinfl R)A*f u ) u , (9) 

^ 3^/2RRR 2 y 1 y ' V ; 

From the non- vanishing Weyl scalars fa, fa and tp4, it is observed that the rotat- 
ing non-stationary de Sitter model (3) is an algebraically special (precisely, type II: 
C a f }C ^£j l ]£ b £ c = 0, with ifio = fa = [9]) in the Petrov classification of space-time with 
a null vector t a (4), which is geodesic, shear free, expanding (6 = \t a a = rR~ 2 ) as well 
as non-zero twist {Co 2 = ^([ a -b]^ a ' b = a? cos 2 Q R~ 2 R~ 2 ) . Here the function A*(u) does 
not involve in the expression of expansion 9 and twist Co 2 of the null vector £ a . This 
means that the physical properties of this null vector l a are same for both stationary 
[2] as well as non-stationary (3) rotating de Sitter models, though they have different 
gravitational fields with different energy momentum tensors. 

The expressions of ip^ and tp4. above involve the derivative of A* (it) as A*(u,) M , 
coupling with the rotational parameter a. So at some point when A* (it) sets to be a 
constant for non-zero rotation (a / 0), both fa and ^4 will vanish. At that point the 
gravitational field of the observer will be of type D (fa ^ 0) in the Petrov classification 
of stationary space-time. That is, the space-time becomes the rotating stationary de 
Sitter solution with cosmological constant A* [2]. 
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The metric (3) can be expressed in the coordinate system (t, x, y, z) as 

ds 2 = dt 2 - dx 2 - dy 2 - dz 2 

{r 6 A*(i r)} r 1 1 1 2 

+ , a , — dt - -r-2- — 2T {r(x<ix + ydy) + a(xdy - ydx)} zdz , 

3(r 4 + a 2 z z ) L [r z + a z ) r J 

where r is defined, in terms of x, y and z [9] 

4 / 2 , 2 , 2 2\ 2 22 n 

r — (x + y +z — a jr — a z = 0, 

with the following transformations 

x = (r cos0 + a sinc/>) sin#, 
y = (r sin0 — a cos0) sin#, 
z = r cos6>, t = u + r. 

Then, the above transformed metric is in the Kerr-Schild form with 

gH = r) ab + 2H(t,x,y,z)£ a £ b , (10) 

where r] a b is the flat metric and 

r 6 A*(t,r) 
■ti(t, x, y, z) — -. j ■ 5 pT-, 
3(r 4 + a z z z ) 

t a dx a = dt — jtt- {r(xdx + ydy) + aixdy — ydx)} zdz. 

[r A + a z ) ' r 

In u-coordinate system the null vector £ a is given in (4) above. The Kerr-Schild 
ansatz (10) confirms that the rotating non-stationary de Sitter model (3) is a solution 
of Einstein's field equations of non-constant curvature. 

The rotating non-stationary de Sitter solution has an apparent singularity when 

A* = r 2 — r 4 A*(u)/3 + a 2 = 0. This equation has four roots r + + , r + _, r_ + and r 

They are found as 



r ±(±) 



2A*( 



1 {3±^9 + 12a^A*(u)}. (11) 



Now let us denote these roots r + + , r + _ , r_ + , r as n, rs, respectively (for 

the simplicity, rj, j = 1, 2, 3, 4). Then, these roots have the following relation 

i \i \i \i \ 3 r 2 r 4 A*(u) 2 1 
(r - n)(r - r 2 )[r - r 3 )(r - r 4 ) = - jr h a j. 

Then each root represents the location of the cosmological horizon for the observer 
and associates an area of the horizon at each point at r = rj, j = 1, 2, 3, 4, 

Aj = / Jgle944> dd d(j) 

Jo Jo r=r i 
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= 4yr{r| + a 2 }. (12) 

According to Bekenstein-Hawking area-entropy formula [3], these areas Aj will de- 
termine the entropies Sj of the horizons of the de Sitter model (3) by the relation 
Sj = Aj/4 [3]. Thus, we find them as 

Sj = ir{r] + a 2 }. (13) 

The gravity of the cosmological horizons is determined by the surface gravity, defined 
by Kn a = n h Vi ) n a in [5], where the null vector n a given in (4) above is parameterized 
by the coordinate u, such that d/du = n a V a , and has the normalization condition 
t a n a = 1 with the null vector £ a . The surface gravities Kj associated at each r = rj 
are found as below: 

k p = ^A*(u){r p + r p+1 }{r p+1 -r p } for p = 1,3, (14) 
% A*(n){r p _! + r p }{r p _! - r p } for p = 2,4, (15) 



m 2 



where R 2 = r| + a 2 cos 2 0, j = 1,2,3,4. The surface gravities Kj may be regarded as 
the gravitational field on the cosmological horizons. From (14) and (15), we observe 
that k\ and K2 will be zero when n and r<i coincide. Similarly, K3 and K4 will vanish 
when r3 = r±. 

The coincidence of two roots n and T2 leads to a condition that {9 + 12a 2 A* (u)} = 
0. This condition implies that r% and also coincide. Then all roots take the form 
r\ = T2 = — r3 = — V4 = A /{3/(2A*(u))}. Accordingly, the area of the horizon at each 
point ri, r2, r3 and are found as 

Aj = (16) 
J A*(u) 

This implies that the entropy associated with each point becomes 

S j = ——, j = 1,2,3,4. (17) 
3 4A*(u) 1 ; 

From this expression we find that the entropies Sj at rj are inversely proportional to 
the cosmological function A*(n). It is to mention that the value of A* (it) is supposed 
to reduce according to the retarded time u change. Consequently the entropies 5^ 
may increase, as the function A*(n) reduces. It is found that once the function A* (it) 
takes the constant value, as in the case of stationary rotating de Sitter universe [2], 
the entropies Sj associated with rj will take constant values. 
The angular velocities f2j for the horizons are found as 

/ Qiid>\ arfA*(u) 
n = lim = \> is 
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The coincidences of the roots (ri = r2, and = r±) imply that rj = 3/(2A*(u)). Then 
the angular velocities take the forms 

An 

= - t a*(«). (19) 

This indicates that angular velocities are directly proportional the cosmological func- 
tion A*(u), and the effect of the change in A*(u) will certainly affect on the angular 
velocities associated with rj. It is also emphasized that the angular velocity Qj given 
in (18) will vanish when the rotational parameter a tends to zero, showing the fact 
that there is no angular velocity for non-rotating de-Sitter space-time with vanishing 
rotational parameter a = 0. 

It is quite interesting to discuss the nature of the non-rotating (a = 0) non- 
stationary de Sitter model with A*(u). Although the non-rotating de Sitter model 
does not explain the complete structure of the space-time, the metric is very simple 
without much mathematical expressions. Thus, when one sets the rotational parame- 
ter a to zero, the metric (3) reduces to non-rotating de Sitter model as 

ds 2 = {l - -r 2 A* («)} du 2 + 2dudr -r 2 (d0 2 + sin 2 'Ode/) 2 ). (20) 

In this situation, the Weyl scalars ^2, ^3 and 1^4 given in (9) are vanished showing that 
the space-time becomes the conformally flat (C abc d = 0). Then the energy-momentum 
tensor (5) takes the form 

KT ab = ~rA*{u), v £ a £ b + A*(u)g ab (21) 

with its trace KT = 4A*(n). Here the energy-momentum tensor (21) involves a 
Vaidya-like null radiation term — ^rA*(u), u £ a £b, which will vanish when r — > 0, and 
satisfies the energy conservation equation T ab . b = 0. However, it still maintains the 
non-stationary behavior A*(u) 7^ constant, showing that the space-time of the observer 
is naturally time dependent even at r — > 0. The energy momentum tensor (21) will 
become the one of the original de Sitter model [1] when A*(u) takes a constant value. 
Using the energy momentum tenor (21) we find Einstein's field equations G ab = —KT ab 
for the non-rotating metric (20) as follows 

Rab~\R 9ab + A{u)g ab = -Ti b NS) , (22) 

where the non-stationary evolution part of the energy-momentum tensor (20) is given 
by 

which has zero-trace T^ NS ^ = 0. It is to emphasize that the universal constant K does 
not involve in the field equations (22). For future use we also present the null energy 
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density p,*, energy density p* and pressure p of the non-rotating de Sitter metric as 
follows 



M * = -— A*(«) )tt p* = 



A*(«) 



P = - 



A*(«) 



(23) 



A' 



From (23) we find the equation of state w = p/p* = — 1 with the negative pressure 
of variable A* (it). This shows the fact that our non-stationary de Sitter solution (20) 
is in agreement with the cosmological constant (A*) de Sitter solution possessing the 
equation of state w = — 1 in the dark energy scenario [10, 11, 12]. 

The metric (20) has an apparent singularity with a horizon at r± = ijSA*^ 1 )^)} 1 / 2 . 
The entropies S± and the surface gravities k± of the horizon associated with r = r± 
for the non-rotating metric (20) are found as 



From these expressions we observe that S± are inversely proportional to A*(u), whereas 
k± are not. It is emphasized that, according to the change of the retarded time u, 
the value of A*(u) may decrease. This ensures that the entropies S± of the non- 
stationary de Sitter model will increase, whereas the surface gravities n± decrease. 
Such observations of changes of the values of S± and n± can be found only in the case 
of non-stationary de Sitter model (20). This is the important aim of the study of non- 
stationary de Sitter space-time. Our results of the non-rotating de Sitter solution are 
in agreement with those of Gibbons and Hawking [3], when the cosmological function 
A*(u) tends to a constant A*. 

The Kretschmann scalar for non-rotating de Sitter model (17) takes the form 



which does not involves any derivative term of A* (u) , and will not change its value at 
r —>■ and r — > oo. The above Kretschmann scalar will become the one of original de 
Sitter model when A*(u) takes a constant value. That is, though the energy momentum 
tensor (21) for non-rotating non-stationary model is found different from the one of 
non-rotating stationary de Sitter solution, the forms of Kretschmann scalar for both 
non-rotating non-stationary (20) and stationary models [1] have similar structures 
with a difference in nature of the cosmological function A*(u). 

Our result discussed here includes the following cosmological models: (i) original de 
Sitter when a = 0, A*(u)= constant [1], (ii) rotating stationary de Sitter when a ^ 0, 
A*(tt)= constant [2], (hi) non-rotating non-stationary de Sitter when a = 0, A*(u) / 
constant, (iv) rotating non-stationary de Sitter when a^0, A*(u) / constant. From 
the study of non-stationary models it is observed that the gravitational field of the 
space-time of the rotating model (3) is algebraically special in the Petrov classification 
whereas the non-rotating one (20) is conformally flat with the energy momentum 




and 




(24) 



y — o riabcd A * 

A" = -Kabcd-K — —A 




(25) 
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tensor describing non-empty space. It is hoped that known stationary works, rotating 
or non-rotating, with cosmological constant, may be extended to the non-stationary 
ones by using the non-stationary de Sitter model, rotating (3) or non-rotating (20). It 
is to note that to the best of the author's knowledge, these non-stationary de Sitter 
cosmological solutions (3) and (20) have not been seen discussed before. 
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Appendix A: Energy-momentum tensor conservation equa- 
tions 

In this appendix we establish the fact that the energy-momentum tensor (5) satisfies 
the conservation equations T ab . b = 0. These are four equations, which can equivalently 
be expressed in three equations - two are real and one complex, using Newman-Penrose 
spin coefficients. Hence, we find the following 

Dp* = (p* + P )(p + p) +UJK* +LUK*, (Al) 

Dp* + Vp* + 5uj + Su = p* {(p + p) - 2(e + e)} - (p* + p)(p + p) 

-u(2TT + 2(3-f)-u(2Tt + 2(3-T), (A2) 

Duj + 5p = p*K* + (p* + p)(r - 7f) + uia - u{2e -2p-p) (A3) 

where k*, t, tt, etc. are spin coefficients, and the derivative operators are defined by 

D = td a , V = n a d a , 5 = m a d a , 5 = fh a d a . (A4) 

These are general equations for an energy- momentum tensor of the type (5). 

Now, in order to verify whether components of T ab with the quantities p*, p*, 
p and co given in (6) for the rotating non-stationary de Sitter metric satisfy these 
conservation equations (A1-A3) or not, we present the NP spin coefficients for the 
metric (3): 

K* = a = X = e = 0, 

1 A* 
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(2ai - R cos 9) cot 9 



2^/2RR sm.9 ' 2^/2i?' 
iasin# iasin# 
7T = ^2RR> T = "72i^' (A5) 
1 
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{r--r 2 A*(u)}R- A* 



^ = — ; — ^ „o zar 4 sin 9A*(u) u 
6^2 RR 2 v J ' 

The derivative operators (A4) are given as follows: 

D = d r , 

^ = jp{(r 2 + a 2 )d u -^-d r + ad <p }, (A6) 

If 2 1 

5 = -rr^ i a sin 9 d u + <9 e + ^— k 

where A* = r 2 — r 4 A*(u)/3 + a 2 . The equations (Al) and (A3) are comparatively 
easier to verify than (A2). Therefore, we shall not show their verification here except 
for the equation (A2). Now, by virtue of (6) and (A6), the left side of (A2) is found 
as 



+ 



r 5 a 2 sin 2 , 2r 3 a 2 cos 2 6>A* 
FT 

{■ 



2r 6 a 2 r 4 



+ 



3KR 2 R 2 R 2 3KR 2 R 2 R 2 
2a 2 r 4 sin 2 (9 



^WWWr 2 ( 3a cos 9 + r ) l A {ulu - (A7) 



which can be shown equal to the right side of (A2), by using (6) and (5A). This leads 
to the conclusion of the verification that the energy-momentum tensor (5) satisfies the 
condition T ah . b = 0. 
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